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Based on analytical and numerical approaches, we investigate thermodynamic properties of 
CuSe 20 s at high magnetic fields which is a candidate for the strong intra-chain interaction in 
quasi one-dimensional (ID) quantum magnets. Magnetic behavior of the system can be described 
by the ID spin-1/2 XXZ model in the presence of the Dzyaloshinskii-Moriya (DM) interaction. Un¬ 
der these circumstances, there is one quantum critical field in this compound. Below the quantum 
critical field the spin chain system is in the gapless Luttinger liquid (LL) regime, whereas above it one 
observes a crossover to the gapped saturation magnetic phase. Indications on the thermodynamic 
curves confirm the occurrence of such a phase transition. The main characteristics of the LL phase 
are gapless and spin-spin correlation functions decay algebraic. The effects of zero-temperature 
quantum phase transition are observed even at rather high temperatures in comparison with the 
counterpart compounds. In addition, we calculate the Wilson ratio in the model. The Wilson ratio 
at a fixed temperature remains almost independent of the field in the LL region. In the vicinity of 
the quantum critical field, the Wilson ratio increases and exhibits anomalous enhancement. 

PACS numbers: 75.10.Jm; 75.10.Pq 


I. INTRODUCTION 

From the theoretical point of view, the behavior of 
quantum magnets usually is determined by the Heisen¬ 
berg model of interaction between spins. Moreover, it is 
known that a lack of symmetry in some quantum magnets 
leads to consider an extra interaction which is known as 
the Dzyaloshinskii-Moriya (DM) interactionii^. In fact, 
the DM interaction arises from the spin-orbit coupling 
and generates many surprising characteristics. 

Experimentally, some quasi-lD antiferromagnetic 
systems are known to be described by the DM 
interaction^rJ^. Recently a compound with the chemical 
formula, CuSe 20 ^, is known as a candidate for the quasi- 
ID quantum magnets with the DM interactioni^P— . DM 
interaction in a realistic model with multiorbitals on the 
ligand oxygen ions is the origin of anisotropic behavior 
which exists in this compoundi^ First, the structure of 
CuSe20^ is characterized by chains with Cu ions aligned 
along the c axis of a monoclinic latticei^. The mentioned 
crystal structure is reinvestigatedi^ in 2006 and to much 
higher precision is confirmed. By investigating the mag¬ 
netic susceptibility, a strong intra-chain antiferromag¬ 
netic exchange interaction is revealedi^. In addition to 
the mentioned intra-chain exchange coupling, it is argued 
that a non-frustrated very weak inter-chain coupling ex¬ 
ists between the nearest-neighbor spinsii. This interac¬ 
tion causes a 3D antiferromagnetic ordering in the range 
of temperature below 17 K. Using the Raman scattering 
study, it is confirmed that CuSe 20 ^ is characterized by 
a moderate, non-frustrated inter-chain couplingi^. Using 
the electron spin resonance method, the low-temperature 
behavior is also studiedi^. It is found that the symmetric 
anisotropic exchange and the antisymmetric DM interac¬ 
tion are almost the same in this system. In a very recent 
work, using the bulk magnetization, neutron diffraction, 
muon spin relaxation and antiferromagnetic resonance 


measurements the magnetic behavior of this compound 
is studied. The long-range Neel order below Tm = VI K 
is characterized^^. 

Theoretically, Using the high-temperature series 
expansion^^ and the Bethe-ansatz approach the tempera¬ 
ture behavior of this compound is studie d^^d^i^° . A very 
good agreement with the experimental results is found in 
the temperature region T > Tjv = 17 AT where the ther¬ 
mal fluctuations are dominant to the effect of the inter¬ 
chain exchange interaction. These theoretical studies are 
limited to a range of the magnetic field less than 5 T. 
Since the intra-chain exchange interaction is very strong 
in this compound (J = 157 K), the saturation field is 
very high and experimentally is not available now. 

Since the physics of the mentioned compound can be 
obtained by the ID spin-1/2 antiferromagnetic Heisen¬ 
berg model with added the DM interaction, the inves¬ 
tigation of the high-held behavior is possible from the 
theoretical point of view. Therefore, in this paper, by 
using the analytical fermionization approach and the nu¬ 
merical stochastic series expansion QMC method (the 
ALPS^ code), we draw a theoretical clear picture of the 
probable behavior of the compound CuSe 20 ^ in very 
high magnetic helds. Specially, we study the Wilson ra¬ 
tio in this compound. Since, recently the Wilson ratio 
on a gapped spin-1/2 two-leg ladder has been measured 
in an experimental, we suggest that the CuSe 205 com¬ 
pound is a very good candidate for measuring the Wilson 
ratio. In addition, this is an opportunity for expanding 
the knowledge on the universal nature of ID quantum 
liquids. 

The paper is organized as follows. In the forthcom¬ 
ing section we introduce the model and map it into an 
effective ID spin-1/2 XXZ model in the presence of a 
magnetic field. Using the fermionization approach we di¬ 
agonalize the Hamiltonian. In section HI, we present our 
analytical and numerical results and compare them. Fi- 
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nally, we conclude and summarize our results in section 

IV. 


II. THE MODEL 


The Hamiltonian of the ID spin-I/2 antiferromagnetic 
Heisenberg model with added DM interaction in a mag¬ 
netic field is written as 

N N 

n = + X ^.+i) 

N 

- gtisHYS], (1) 

i=i 

where Sj is the spin-I/2 operator on the j-th site, J > 0 
denotes the exchange coupling constant, H is the applied 
magnetic field and = Dz is known as the DM vector. 
The exchange coupling and the DM vector are measured 
J = 157 K, D = 0.05 J respectively in the CuSe 205 
compound. At first, by considering a rotation^^ around 
the Z axis as exp(=Fio) in which tana = 

the Hamiltonian is transformed to the ID anisotropic 
spin-1/2 XXZ model in a longitudinal magnetic field 


u = jY.(sjsu, + s]sj„ + jS;su,) 

N 

- gfisHYS], (2) 

with an effective exchange interaction J = \/+ D"^. 
Since the exchange anisotropy is less than the XY 
anisotropy (J < J), the ground state of the system is 
in the Luttinger liquid (LL) phase at zero temperature^!. 
By increasing the magnetic field from zero, up to the crit¬ 
ical field He = 2+2. ground state remains in the LL 
phase where a quantum phase transition into the ferro¬ 
magnetic phase with saturation magnetization along the 
field will occur. 

Theoretically, the energy spectrum is needed to study 
the thermodynamic behavior of the model. An efficient 
method for diagonalizing is the analytical fermioniza- 
tion technique. At the second step, by implementing the 
Jordan-Wigner transformations. 


^3 ~ 2 ’ 

Y = a]exp(i7r^aj’ai), 

l<3 

S~ = aj exp(-i7ry^ aj’g;), 

l<3 


the transformed Hamiltonian is mapped into a ID model 
of interacting spinless fermions, 




NgueH NJ 
2 ~ 




2 ^ 


{a]aj+i + a^aj) 


j 

(J -f gg,BH) ■ ( 4 ) 
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Then, using the Wick’s theorem, the fermion interaction 
term is decomposed by some order parameters which are 
related to the spin-spin correlation functions as^^ 

71 = {a]aj), 

72 = {a]aj+i), 

T's = (OiSVi)' (5) 

By utilizing these order parameters and performing a 

Fourier transformation as aj = 

also using the following unitary transformation 

ttfe = cos{k)l3k - isin{k)/3l, (6) 

the diagonalized Hamiltonian is given by 

^/ = E - i). (7) 

k— — T^ 


Where the dispersion relation is 

e{k) = \/a{ky + Ab{k)^, 

a{k) = {J — 2^2 J)cos{k) + 2ji J — {gjisH + J), 
b{k) = Jj3sin{k). ( 8 ) 

In order to solve TLf, the following equations should be 
satisfied self-consistently 

^ 1 r a{k), 1 1 , 

e(fc)^l + e"W/^BT 2 ^’ 

1 r., 1 ^ 

T'2 2Tr + 2 ^''°®^^^’ 

1 r „ b(k), 1 lx 

^ ^(1 + eeW/i.BT - (9) 

For this purpose, we have written a computational code 
to solve the self consistent equations. Using our code, 
first we find a set of 7 i, i = 1, 2, 3 with fixing the values 
of the exchange coupling, the DM vector and the mag¬ 
netic field. Then, by including 7 ^,* = 1,2,3 into the 
dispersion relation, the energy spectrum of the model is 
obtained. Finally the thermodynamic functions such as 
the free energy, the specific heat, the magnetization and 
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the susceptibility are obtained as 
KbT 


f 

Cv 


27r 


dk ln(2cosh{—-^-\—)), 


1 


m = 


X = 


SttKbT^ 

27r J 

dm 




'2KbT' 
dk (e{kfsech^{^^)), 
^ _ a{k) 


e(fc)^l + e®(fc)/i^flT^ 2e(fc)'’ 


( 10 ) 


The analytical fermionization approach is used for in¬ 
vestigating the behavior of the system from a theoretical 
point of view. 


III. RESULTS 

In this section we present our results regarding the 
thermodynamic behavior of the ID spin-1/2 XXZ model 
in the presence of a longitudinal magnetic field in with 
DM interaction. As mentioned before, the field of our 
study is the CuSe 205 compound at high magnetic fields. 
To this end, we use the experimental data in our calcu¬ 
lation, for compatibility between theory and experiment. 
The maximum size of spin chains in QMC simulations is 
100 particles. The QMC simulation is performed for spin 
chains under the periodic boundary condition with the 
maximum 3x10® equilibration sweeps and 6x10® mea¬ 
surement steps. It is believed that the simulation ap¬ 
proaches are very good to find the results in cases where 
the performing of a real experiment is impossible. Since, 
there are not any real experimental results for the behav¬ 
ior of the CuSe 20 ^ compound in high magnetic field, we 
have used the QMC simulation. We have also solved 
the integral equations (Eq. ([5])) numerically to find the 
energy spectrum of the model in the fermionization ap¬ 
proach. Figs.[lja) and 1(b) show the magnetization ver¬ 
sus the longitudinal magnetic field at different values of 
the temperature by using fermionization approach and 
SSE QMC simulations respectively. At very low temper¬ 
ature, both of them show that the magnetization starts 
from zero, and it confirms that the excitation spectrum 
of spin chain is in the gapless LL region. By increas¬ 
ing the magnetic field, the magnetization increases very 
slowly because of the strong exchange interaction. Thus, 
the very large value of the field He = 2300 KOe has the 
role of quantum critical field and in the region of H > He 
the magnetization will be saturated. The spin chain sys¬ 
tem undergoes a crossover from the gapless LL phase to 
the gapped saturation magnetic phase. As is seen, a very 
good agreement exists between the QMC simulation and 
fermionization approach results. By increasing the tem¬ 
perature, the saturation plateau will be destroyed due to 
the thermal fluctuations. Experimentally the magnetiza¬ 
tion process was studied in the temperature T = 4.2 K 
and the magnetic fields up to SOATOe^. Fig.[2]shows our 


analytical and numerical results together with the exper¬ 
imental datap2£. It is clearly seen that, at T = 4.2 K our 
theoretical results have a very good agreement with the 
experimental data. 

Figs. [3] (a) and 3(b) illustrate the molar susceptibil¬ 
ity versus the magnetic field with both SSE QMC and 
fermionization approaches. At quite low temperature 
(T = l.SAT), by increasing the magnetic field up to 
2300 KOe a sharp peak emerges in x(iL), which is an 
indication of the quantum critical field. Due to the very 
large exchange interaction, as depicted in Figs. [3] (a) and 
3(b), at magnetic fields nearly up to 2000 KOe the molar 
susceptibility is the same within all the considered range 
of temperature which is one of the important signatures 
of the LL phase. By more increasing of the magnetic 
field up to the quantum critical field in the region of LL 
phase the molar susceptibility decreases with increasing 
temperature. As expected, enhancing of the temperature 
causes the sharp peak to fade. In other words, the effects 
of quantum fluctuations disappear and thermal fluctua¬ 
tions are dominant. The behavior of molar susceptibility 
in the gapped saturation magnetic phase is different from 
the gapless LL phase. In other words, we observe that the 
molar susceptibility increases with enhancing of temper¬ 
ature. To complete our study, we calculate the specific 
heat for different values of the temperature and the mag¬ 
netic field. As usual, figures have been plotted with two 
numerical and analytical results. Figs. |4Da)-4(d) show 
the magnetic specific heat versus the temperature. Be¬ 
low the critical field H < He, with lowering the temper¬ 
ature, there is a broad peak about 100 K which is called 
a Schottky peak, then with more decreasing temperature 
we observe a second peak in the specific heat and comes 
down to zero linearly, indicating the existence of the LL 
gapless regime. By increasing the magnetic field and near 
to the quantum critical field the height of the second 
peak decreases and then vanishes. Figs. S] (c) and 4(d) 
have been plotted for H > He- Both QMC simulation 
and fermionization approaches give us the same results. 
Above the quantum critical field H > He = 2300 KOe, 
a crossover occurs into the gapped saturation magnetic 
phase and there is no trace of the second peak. On the 
other hand, Cv{T) remains zero up to a threshold tem¬ 
perature, which is known as an indication of the energy 
gap in the saturation magnetic field. Next, we consider 
the specific heat versus the magnetic field C'y(iL) at a 
fixed temperature. As illustrated in Figs. [5] (a) and 5(b) 
the curve of Cy (H) has been plotted at various fixed tem¬ 
peratures T = 20 K, 30 K, 40 K and 50 K. It is clearly 
seen from this figure that the specific heat at low tem¬ 
perature, shows two-peaks structure on both sides of the 
quantum critical point in good agreement with results 
obtained for the Kagome antiferromagnet chains2^ and 
Ising model^S. Within the temperature mentioned above, 
there is a minimum at He = 2300 KOe, which is an indi¬ 
cation of the existence of the zero-temperature quantum 
phase transition. This minimum fades with enhancing 
temperature, because of thermal fluctuations. The sig- 
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nature of existence of zero-temperature quantum phase 
transition in the curve of Cv{H) even in this range of 
temperature (up to 100 K) is because of the strong ex¬ 
change interaction in this compound. Based on these 
results, it seems that quantum correlations^^ as the en¬ 
tanglement in the CuSe 205 compound can be detected 
experimentally at high temperatures which can be con¬ 
sidered as an open problem. In addition, we have fo¬ 
cused on the theoretical evaluation of the Wilson ratio 
in the CuSe20^ compound. It is known that the Wilson 
ratio^. 


Rw 


3 ^ 5MB a/T’ 


( 11 ) 


is a crucial parameter for characterizing the LL region. 
In this relation, y is the magnetic susceptibility and Cy 
is the specific heat. In principle, the Wilson ratio quanti¬ 
fies the interaction effect and the spin fluctuations which 
enhance the magnetic susceptibility. It was illustrated 
that Rw is 2 for the spin-1/2 Kondo lattice in the single¬ 
impurity limit^. It will be 1 for a system of noninteract¬ 
ing electrons in a metal^ and independent of the tem¬ 
perature. Since the spin fluctuations are enhanced in 
strongly correlated systems, Rw > 1 is expected . Since 
the Wilson ratio has recently been measured in exper¬ 
iments on a gapped spin-1/2 Heisenberg ladder—, the 
study of the Wilson ratio in the ID spin systems has also 
attracted much interest. At zero temperature, Rw ex¬ 
actly equals 2 for an isotropic antiferromagnetic spin-1/2 
chains^. In this system, Rw also depends on the tem¬ 
perature. Recently the effect of the polarization on the 
Wilson ratio of Fermi gases in ID has been studied^^— . 
It is found that Rw exhibits anomalous enhancement at 
critical points. 

Here we have calculated the Wilson ratio in our model. 
In Figs. [6] (a) and 6(b), Rw has been plotted as a 
function of the magnetic field at a fixed temperature 
T = 20AT. The special value of 20AT has been selected 
since CuSe20z compound behaves as a quasi-lD system 
at temperatures higher than Tjv = IVAT. Our analyti¬ 
cal and numerical results show that the Wilson ratio is 
2 < Rw < 3 and independent of the magnetic field in 
the region oi H < He, where the LL phase is expected. 
An anomalous enhancement of the Wilson ratio in the 
vicinity of the quantum critical point is seen. This phe¬ 
nomenon has also been observed near the metal-insulator 
transition in simulations of a three dimensional quantum 
spin liquid^. Here, this anomalous divergence is mainly 
an indication of the quantum critical point. In addition, 
this behavior confirms that the susceptibility of the sys¬ 
tem increases in the vicinity of the quantum critical point 
much faster than the specific heat. In other words, the 
critical exponent of the susceptibility should be larger 
than the critical exponent of the specific heat. It is inter¬ 
esting, what happens when we change the value of DM 
interaction, however the value of measured DM vector is 
constant {D = 0.05J). Since the transformed Hamilto¬ 
nian (Eq. (2)) is a ID spin-1/2 XXZ model, we expect 


that by increasing the DM interaction the model remains 
in the LL phase and no quantum phase transition occurs 
at zero temperature. Among the thermodynamic func¬ 
tions, the evaluation of the behavior of the Wilson ratio 
in respect to the DM vector is interesting. In Figs. [7] , 
the Wilson ratio has been plotted as a function of the 
DM interaction at a fixed temperature, T = 20K and H 
= 500KOe. As is seen, Rw remains almost constant and 
independent of the DM vector in this region. We expect 
this behavior, since the system remains in the LL phase 
by increasing the DM vector. 

Finally, using the mean field approximation^^ we 
have estimated the critical Neel temperature in the 
CuSe20^ compound. We have approximated the 
CuSe 20 ^ compound with a system of coupled parallel 
XXZ chains where the exchange interaction among neigh¬ 
boring chains is Jm = O.IJ = 15.TA'. The condition for 
determining Tjy is where z is the co¬ 

ordination number of the lattice. Using this relation, we 
found the Neel temperature Tjq ~ 22 AT, which is in good 
agreement with the experiment. 


IV. CONCLUSION 

In summary, we have investigated the high magnetic 
field behavior of the quasi-ID CuSe20z compound. It 
has been suggested that the behavior of this compound 
can be obtained by the ID spin-1/2 Heisenberg model 
with the Dzyaloshinskii-Moriya (DM) interaction. At the 
first step, we have mapped the Hamiltonian of the system 
into an effective ID spin-1/2 XXZ model. Then, analyt¬ 
ically the effective Hamiltonian is diagonalized using the 
fermionization approach. Using the energy spectrum, the 
thermodynamic functions are calculated. In addition, the 
numerical SSE QMC simulation is also applied and the 
same thermodynamic functions are calculated for finite 
chain systems up to V = 100 spins. 

Results show that the system is in the LL phase in the 
presence of the external magnetic fields up to 2300ArOe. 
With further increase of the magnetic field, the system 
undergoes a phase transition into the saturation mag¬ 
netic phase. The molar susceptibility is independent of 
the temperature in the LL phase. Also, specific heat ver¬ 
sus temperature shows two-peaks in the LL region, which 
is a very important signature of the LL phase. In addi¬ 
tion, in the figure of the specific heat versus magnetic 
field, two-peaks structure on both sides of the quantum 
critical point can be observed in high values of tempera¬ 
ture ~ 50Ar. The mentioned behavior suggests that the 
quantum correlations in this compound should be very 
strong and their effects can be observed at high temper¬ 
ature region. In the topic of ” quantum entanglement”— 
and its applications, having the quantum correlations at 
high temperature region is a challenge. Our results show 
that the CuSe 20 ^ compound can be a very good candi¬ 
date for studying the entanglement in the experiment. 

We also tried to find an estimation of the Wilson ratio’s 
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order in the CuSe20^ compound. Our results show that 
the Wilson ratio should be in the range 2 < < 3. On 

the other hand, it remains almost constant up to very 
high magnetic fields ~ 2000KOe. 
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FIG. 1: Magnetization versus magnetic field at zero temper¬ 
ature showing a critical field at He = 2300 KOe (a) fermion- 
ization approach (b) SSE QMC simulation. 
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FIG. 2: Magnetization versus magnetic field at temperature 
T = 4.2 K which is coincident to experimental data. Fermion- 
ization approach and SSE QMC simulation has been carried 
out. 




FIG. 3: Molar susceptibility versus magnetic field at different 
temperatures. The asymptotical behavior in the vicinity of 
zero temperature is an indication of cross-over to a saturation 
magnetic phase (a) fermionization approach (b) SSE QMC 
simulation. 
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FIG. 4: Magnetic specific heat versus temperature below 
critical magnetic field = 2300 KOe (a) fermionization 
approach (b) SSE QMC simulation, (c) Magnetic specific 
heat versus temperature above critical magnetic field He = 
2300 KOe with fermionization approach (d) SSE QMC sim¬ 
ulation. 




FIG. 5: Magnetic specific heat versus magnetic field at differ¬ 
ent temperatures (a) fermionization approach (b) SSE QMC 
simulation. The quantum critical field shows its effect at 
rather high temperatures. 



FIG. 6: The Wilson ratio as a function of the magnetic field 
at fixed temperature T = 207F(a) fermionization approach 
(b) SSE QMC simulation. 
























FIG. 7: The Wilson ratio as a function of the DM vector at 
fixed temperature T = 20K and magnetic field H = 500KOe 
(a) fermionization approach (b) SSE QMC simulation. 





